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B Physics with NRQCD: A Quenched Study 

Presented by A. Ali Khan ^* 

^Department of Physics and Astronomy, University of Glasgow, UKQCD Collaboration 

We present results on the spectrum of B mesons and heavy baryons, using a non-relativistic formulation for 
the heavy and a clover action for the light quark. We also discuss B meson decay constants and their dependency 
upon the heavy meson mass. 
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1. INTRODUCTION 

For some years, particle theorists have put con- 
siderable effort into calculations of the B meson 
decay constant fs- This quantity is required, to- 
gether with the bag parameter Bb, for a deter- 
mination of the CP-violating phase in the CKM 
matrix. However, relativistic quark methods need 
some reinterpretation if they are to work for h 
quarks [1]. Simulations are generally done at 
lighter quark masses and results extrapolated to 
the h, assuming an expansion in the inverse heavy 
quark mass, as predicted by Heavy Quark Ef- 
fective Theory (HQET), up to quadratic terms. 
Some difficulties in matching with the static re- 
sults are often encountered. Extrapolation has 
also to be performed to obtain the spectrum 
around the B, and it is hard to get sensible val- 
ues for the hyperfine splitting (a successful static 
calculation is in [2]). However, with NRQCD it is 
possible to simulate in mass regions around and 
above the h quark. NRQCD is an effective theory 
which systematically includes higher dimensional 
operators in an expansion in the inverse heavy 
quark mass I/itIq. One can study contributions 
of operators occuring in HQET separately. We 
are performing a calculation of the spectrum and 
decay constants of heavy-light mesons. Our study 
of the spectrum comprises S and P wave states, 
as well as the Aj baryon. P states have recently 
been discovered experimentally and provide a way 
to determine the lattice spacing from the B sys- 
tem itself. 



2. THE COMPUTATION 
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The NRQCD Lagrangian in this calculation 
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and in [3] is corrected through 0(1/171%) 
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where Q is the Pauli spinor describing the heavy 
quark. The gauge links are tadpole improved: 



(2) 



so we use the tree level coefficients for the I/itIq 
terms in (1). For the calculation of the heavy- 
light decay constants it is important to also 
include the l/nig corrections to the currents. 
These can be obtained by doing an inverse Foldy- 
Wouthuysen transformation on the heavy quark 
spinor: 
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Thus we calculate the decay constants from ma- 
trix elements (0|C|M), where 



C = qhTQ 
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is given by the time component of the axial vector 
current (M = B, T = 7570) or the vector current 
(M = 5*, r = 7^), respectively. 

We simulate heavy quarks with masses = 
1.71, 2.0, 4.0, 8.0 and in the static approxima- 
tion. We use light quarks with a tadpole im- 
proved clover action at k values 0.1370 and 
0.1381. The calculation is being done on 
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Table 1 

Bare ground state energies (heavier k). 



aAEp 

aniQ c = 1 c = 1/u'q 

1.71 0.501(9) 0.516(5) 
2.0 0.506(6) 0.520(4) 
2.5 0.510(6) 
4.0 0.511(7) 0.532(7) 
8.0 0.537(7) 
static 0.524(6) 0.542(8) 



3. ANALYSIS AND RESULTS 

3.1. B mass and decay constants 

To extract ground state energies (£"0) and am- 
plitudes we fit smeared- local (C'sl) and smeared- 
smeared (C'ss) correlation functions simultane- 
ously to a single exponential. For sufficiently 
large times one has: 

C'ss ^ ^|e-"^°*, (5) 



CsL ^ ZsZlc-''^"'. (6) 



quenched configurations, generated by UKQCD, 
at /3 = 6.0, fixed to Coulomb gauge, on a lattice 
volume of 16^ x 48. The results presented here 
are preliminary, based on an ensemble of ca. 50 
configurations at itIq = 8.0 and ca. 60 configura- 
tions at the other heavy quark masses. Using 130 
configurations the lattice spacing has in the calcu- 
lation of [4] been determined from light spectro- 
spcopy (mp) to be = 1.891(103) GeV and Kg 
to be 0.139245(33). These data are compared to 
a previous calculation with clover light fermions 
with clover coefficient c = 1 (i.e. not tadpole 
improved) at k = 0.1432 and 0.1440 and heavy 
masses 1.71, 2.0, 2.5, 4.0 at the same /3 value and 
lattice volume, with a lattice spacing from light 
spectroscopy of 2.05(10) GeV. In the quenched 
approximation experience shows that a depends 
upon the system used to determine it, so it is 
desirable to use a value from the B system it- 
self, like the S — P splitting. Until sufficiently 
accurate simulation results for this are available, 
we use the values from light spectroscopy for our 
heavy-light results. We expect the light quark 
momentum to set the energy scale in B mesons. 
The K values are tuned such that the heavier and 
the lighter ones of each give approximately the 
same tt mass. Implemented are ^Sq, ^Si, ^Po, 
^Pi, ^P'j and ^ Pi states, taking also cross corre- 
lations between P states into account. Smearing 
functions are hydrogen-like ground and 1st ex- 
cited state wave functions. For each meson we 
calculate all combinations of smeared and local 
operators at source and sink. The Aj baryon is 
smeared at the source and local at the sink. 
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Figure 1. Ground state energies for tadpole im- 
proved clover fermions (squares, k = 0.137) and 
clover fermions with c = 1 (bursts, k = 0.1432). 



Having extracted , we can write down the ma- 
trix elements. For the pseudoscalar meson we 
have \/2Zl = a^^^fp\/Mp and for the vector me- 
son ^/2Zl = a^^^ fvV^v ■ The meson mass we 
obtain from the relation 



M = Eo + A 



(7) 



A includes the renormalized heavy quark mass 
and the divergent zero point energy of the quark. 
The values used in this work are calculated by 
C. Morningstar [5]. Ground state ener- 

gies for the heavier k values are listed in table 1. 
We find that the results are larger for the tad- 
pole improved than for the non-tadpole improved 
light fermions, as shown in figure 1. This is ex- 
pected from a perturbative analysis of the effect 
of increasing the coefficient of the clover term. 
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Table 2 



Meson masses and decay constants (tadpole im- 
proved light fermions, k = 0.1370), all values in 
lattice units. 





Mp 




fvVW 


spin avg. 


1.71 


2.25 


0.177(5) 


0.200(5) 


0.194(4) 


2.0 


2.54 


0.186(5) 


0.208(7) 


0.188(5) 


4.0 


4.60 


0.220(7) 


0.232(8) 


0.229(78) 


8.0 


8.17 


0.253(8) 


0.260(10) 


0.258(9) 


static 


CXD 


0.296(9) 
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Figure 2. Lattice matrix element a^l'^ f\/M with 
NRQCD and tadpole-improved light fermions 
(circles: pseudoscalar, squares: vector, diamonds: 
spin averaged, cross: static) plotted against the 
inverse spin averaged meson mass. The lines are 
correlated fits to the NRQCD points. 

Table 3 

Meson masses and fp in GeV. 



c = 1 c = I/mq 





Mp 


fp 


Mp 


fp 


1.71 


4.5(5) 


0.16(4) 


4.1(5) 


0.16(3) 


2.0 


5.1(6) 


0.16(4) 


4.7(6) 


0.16(3) 


2.5 


6.0(8) 


0.18(3) 






4.0 


9.3(7) 


0.19(5) 


8.6(11) 


0.18(4) 


8.0 






15.4(13) 


0.21(3) 


static 




0.28(5) 




0.22(6) 



Thinking of the light quarks non-relativistically, 
the clover term adjusts the effective Darwin term 
in the action and this shifts energies upwards [6] . 
The B/ p ratio is lower for the tadpole improved 
fermions. 

0.3 I \ \ \ \ \ \ \ \ \ \ \ \ r- 
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Figure 3. Chirally extrapolated pseudoscalar ma- 
trix element Zau^^^ fp\/Mp with tadpole im- 
proved light fermions (circles: NRQCD, fancy 
cross: static) and light fermions with c = 1 
(squares: NRQCD, cross: static) and correlated 
fits to the NRQCD points. 

The lattice results for the pseudoscalar, vec- 
tor and spin averaged decay constants are given 
in table 2. The extrapolation of the NRQCD 
points is compatible with the static simulation 
result, as seen in figure 2. This remains true if 
one includes the matching renormalization con- 
stants from the lattice to the continuum effective 
theory. The slope increases for large masses. It is 
considerably higher than previous lattice calcula- 
tions and of the order of QCD sum rule predic- 
tions [7]. This, along with an estimate of the con- 
tributions of the various operators at ©(l/m^), 
is discussed in more detail in [3]. To convert the 
decay constants into physical units, we are calcu- 
lating perturbatively the matching renormaliza- 
tion constants to the full theory. A comparison of 
the chirally extrapolated pseudoscalar matrix el- 
ements Za fp V Mp for c = 1 and c = I/mq clover 
light fermions is shown in figure 3. Chirally ex- 
trapolated values for fp, together with the corre- 
sponding meson masses, are given in table 3. For 
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the B meson we take 0.16(4) GeV as our estimate 
for the quenched decay constant. The major con- 
tribution to this error comes from the uncertainty 
in a. 

3.2. Mass splittings 

The mass splitting AE between hadrons is ex- 
tracted by fitting the ratio of the corresponding 
correlation functions to a single exponential: 

(8) 



SL^^-aAEt 



(0) 
SL 



where C^g-^ denotes the ground state smeared- 
local correlation function of the state and 
^SL ^'''1 correlation function (for the hyper- 
fine splitting) or the Aj correlation function (for 
the Aj— B splitting). 

3.2.1. B* - B splitting 

The fit results for the B* — B splitting at k = 
0.1370 in lattice and in physical units are given 
in table 4. For the lighter k the statistical er- 
rors are considerably larger; there is no visible k 
dependence. So we approximate the splitting by 
our results at k = 0.1370. From HQET we expect 
the splitting to have the following dependence on 
the heavy meson mass: 



AE (X l/Mp. 



(9) 



Performing a correlated fit of the simulation re- 
sults to a linear function, we find that their be- 
haviour is compatible with this expectation. The 
splittings in the mass region of the B are slightly 
lower than the experimental value of 0.46(1) 
MeV, which might be an effect of quenching. 

3.2.2. Ab - B splitting 

The chirally extrapolated results for the Aj — 5 
are listed in table 4. Our splitting is compatible 
with the experimental value 0.36(5) GeV. The 
central value is slightly high, but this might be 
due to contamination with excited states, since 
correlation functions get large errors for larger 
times. 

4. CONCLUSIONS 

This is a report on an ongoing computation of 
the B meson spectrum and decay constants with 



Table 4 

Hyperfine and baryon-meson splittings (clover 
term c = I/mq). 



AE{B* - B) 



AE(At - B) 



m^Q aAE 



AE[MeY] aAE AE[GeY] 



1.71 0.020(1) 

2.0 0.017(1) 

4.0 0.0091(8) 

8.0 0.0042(7) 



38(6) 0.25(7) 0.47(15) 

32(6) 0.25(7) 0.48(15) 

17(3) 0.27(7) 0.50(16) 

8(2) 0.26(7) 0.50(16) 



NRQCD in the quenched approximation. We 
compare results using tadpole improved clover 
light fermions with c = 1 clover light fermions. 
The tadpole improved results turn out to be sub- 
ject to stronger statistical fiuctuations. One ob- 
serves that the tadpole improved ground state en- 
ergies are slightly enhanced. Now also a very high 
quark mass is included into the simulation, and 
we find that the extrapolation of the NRQCD de- 
cay constants to infinite mass is compatible with 
the static value. The matching renormalization 
constants to the full theory are being calculated 
in NRQCD perturbation theory. Using prelim- 
inary perturbative results we quote a value of 
fg = 0.16(4) GeV for the B decay constant. 
NRQCD also provides us with a tool for com- 
puting B* — B and Ai, — B splittings at the B. 
The hyperfine splitting is proportional to l/Mp 
as predicted by HQET. In the region of the B 
mass, we find a splitting slightly lower than ex- 
periment. The Aj — 5 splitting is compatible with 
the experimental value. 

This work was supported by SHEFC, PPARC, 
the U.S. DOE and the NATO under grant number 
CRG 941259. 
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